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1 Introduction 



1.1 Aim of the paper 

In [T| a program was started to create a reference-book on matrix-model r-functions - the new gener- 
ation of special functions, which are going to play an important role in string theory calculations. The 
goal is to extract and considerably extend spectacular results obtained during the golden era of matrix 
model studies in late 80's and during sporadic moments of emerging new interest afterwards (see, e.g., 
reviews [21 [3] and references therein). In a number of steps was made towards realization of this 

program for the most fundamental partition function of Hermitian matrix model. Integrable aspects 
of that theory were earlier considered in ^i3j. Additional progress is made in the QFT-like approach, 
which is being developed in a complementary series of papers [8]. While we are still far from having 
a concise and exhaustive presentation in case of the Hermitian matrix model, time is also coming 
to extend analysis to other matrix-model r-functions. The present paper being a sequel of [1] is a 
direct counterpart of [1] for the Kontsevich [Oj and Generalized Kontsevich Model (GKM) [101 [TT] 
r-functions. For parallel consideration see [12j . Again, integrable and QFT aspects of the theory 
are mainly not included: we concentrate mostly on Virasoro-like constraints, perturbative resolvents 
(multi-densities) and duality properties, which are still insufficiently represented in the literature. In 
fact, the situation with the Kontsevich r- function per se is a little better than with the Hermitian 
model, because its direct relation to the topological field theory [13] stimulated a relatively systematic 
consideration in the past, see [14j-[16j and references therein. Still, the most interesting part of the 
story - that of {p,q)-{q,p) duality - remained almost untouched after preliminary papers [T7 t [T8 l [T9]. 

As soon as the topic can be hardly exhausted within one paper, here we concentrate only on a few 
basic examples leaving further developments for future publications. In particular, we mostly focus 
on calculating resolvents (=loop operator averages) in the Kontsevich model, with a special emphasis 
on its simplest (Gaussian) phase, where exists a surprising integral formula [16], and the expressions 
for the resolvents in genera zero and one are especially simple (in particular, we generalize the known 
genus zero result |20j to genus one, see also |21j). Thus, the program of [T] is realized for the Kontsevich 
model with several important simplifications: in variance with [T\, there are very simple formulas for 
the n-point resolvent for lower genera; there is an integral formula for the Laplace transform of the 
n-point resolvent (summed over all genera!); and the (integral) equation for the generating function of 
all resolvents looks especially simple. In the paper, we give a review of these simplifications in sect. 3, 
while before that, in sect. 2, we discuss various features of generic phases of the Kontsevich model. 
Note that in a generic phase the main method to calculate resolvents is to solve the loop equations 
[22 1 123 ^ [24]. which generally have many solutions, the Gaussian Kontsevich model being the only one 
that has the unique solution. However, there is a counterpart of this unambiguous Gaussian solution in 
the generic case, the solution called Dijkgraaf-Vafa (DV) solution [25]- [27]. Among other features, this 
solution has specific integrable properties. We discuss the DV solution in sect. 4. The generalization 
of results to the Generalized Kontsevich model is contained in sect. 5-6, where we also discuss the p-q 
duality in terms of resolvents and corresponding Riemann surfaces in the example of dualities between 
(2,3) and (3,2) models. 

1.2 Correlators in matrix models 

By essence, the main problem of matrix models one may address to is constructing a quantum field 
theory (QFT) presentation of matrix models. Solving this problem would allow one to effectively deal 
with all other problems. 

The main purpose of QFT study of any model is to evaluate arbitrary correlation functions in an 
arbitrary phase and, after that, to study possible relations ("dualities") between these correlators in 
different phases. In the context of matrix models certain subsets of correlation functions are naturally 
collected into generating functions which will be called resolvents (or multi-densities). They possess, 
at least, three different representations. 

La raison d'etre for (multi-)resolvents is a transparent group-theoretical structure of the Schwinger- 
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Dyson equations (which is obscure in a generic QFT but is immediately obvious in the simple matrix 
models): these are W- (Virasoro in the simplest cases) constraints with a loop-algebra structure. Ac- 
cordingly, correlation functions satisfy the loop-equations, and the loop parameter becomes a natural 
expansion parameter of the generating functions. 

It is still difficult to solve the genuine loop equation and obtain the full generating function, but 
an additional "genus expansion" converts the loop equation into a chain of simpler loop equations 
for partial generating functions, multi-resolvents which can be evaluated straightforwardly one after 
another. Ambiguities arising in this recursive process lead to different sets of multi-resolvents and 
are interpreted as associated with different phases of the theory. This will be our first approach to 
multi-resolvents. 

Multi-resolvents emerge as non-trivial functions of the loop parameter z with singularities of various 
types, both poles and branchings. The second approach deals with them as poly-differentials on an 
auxiliary Riemann surface Sq ("spectral" complex curve), and different phases correspond to different 
choices of the spectral surface and, in addition, to different conditions on the periods of multi-resolvent 
poly-differentials (e.g., if all the periods but the periods of the first multi-resolvent are vanishing, one 
gets to the so called Dijkgraaf-Vafa phase). As usual, poly-differentials on Riemann surfaces are most 
immediately represented as correlators of free fields, hence, this approach is often called conformal 
field theory (CFT) representation. An adequate reformulation of the loop equations suitable for the 
CFT representation is partly worked out in [8l [6] . 

The third approach represents correlation functions via (functional or matrix) integrals. The 
problem, however, is that the multi-resolvents are generating functions of matrix model correlators, 
i.e. derivatives of the matrix model partition function, and, therefore, their representation by integral 
formulas is not a priori obvious. In the Hermitian model this representation is rather straightforward: 
the spectral (or loop) parameter z is introduced through the average of the loop operator 

^ az \ - az ^ ,1, \ - az a 

where (p is the Hermitian matrix that is integrated over. However, a counterpart representation for 
the Kontsevich model remains unclear (see |18j for a very tedious approach to evaluating a few first 
^ for the generalized Kontsevich integral). Worse than that, even integral formulas for the partition 
function are also unknown for most of non-trivial phases of the Kontsevich model. Still, spectacular 
results for correlators of the Gaussian Kontsevich model, due to [16] (discussed in sect. 3), imply that 
the third approach should also be fruitful. Somewhat surprisingly the integral formulas in the Gaussian 
case are most simple not for the multi-resolvents themselves (i.e. not for the quantities subjected to 
the loop equations), but for their Laplace transforms. 

1.3 IZK integral 

Now we specify our general discussion to the case of the Kontsevich model, the main object of the 
present paper. 

The story about the Kontsevich model begins from the Itzykson-Zuber-Kontsevich (IZK) integral 
over n X n Hermitian matrices X 

HMV) = I dX exp (tr AX - tr V{X)) , (1.2) 
depending on the choice of the potential V{x), 

oo 

V{x) = ^SkX^ (1.3) 

fc=0 

and on the background matrix-valued field A. AA(A) here is a normalization factor. This matrix 
model is actually of the eigenvalue type [71 13] : as was first demonstrated by Itzykson and Zuber 
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the integration over angular variables U in X = X^iagU can be done explicitly, leaving the n-fold 
integral over eigenvalues Xi of X in X^iag = diag(xi, • • • , Xn), 

mV) ~ n / '^^^^"''^'''^XTTY . -dft ^""""^ (1-4) 

, J ZA(Aj ij=l,...,n 

1=1 

where \j are eigenvalues of A and A denotes the Van-der- Monde determinant, 

n 

A(x) = ..det xr' = nte-^j) 

Transition from (jl.2|) to (jl.4p is typical for the Harish-Chandra-style character calculus in group 
theory [221 EH]- A particular expansion of the particular IZK integral with pure cubic potential 
V{x) = was related by M.Kontsevich [9] to cohomologies of the moduli space of Riemann surfaces 
and, finally [M], to partition function of topological gravity [13]. Moreover, this particular cubic 
potential case turns out to be related to more general Hodge integrals over the moduli space that 
include A-classes [32j, these latter being related to the Hurwitz numbers [33j. Many properties of the 
integral are, however, independent on particular choice of V{x) and can be addressed in the theory of 
Generalized Kontsevich Model (GKM) p^ [3i]. 

• The first split between different directions of study of the GKM concerns the type of A-dependence 
in (II. 2p . One option is to consider tr AX as a perturbation and represent /(A) as a series in powers of 
tr A'^ with A; > - this is the character phase of the model \30\ |3l] . Instead one can expand around 
a classical solution X = L to the equation of motion V'{L) = A of the full action, then the expansion 
will be in powers of tj. = ^tiL~^ with k > provided the normalization factor J\f in p.2p is chosen 
equal to the quasiclassical value of the integral - this is the Kontsevich phase of the same model. In 
the GKM with monomial potential 

Vp{x) = — - 
p+l 

L is just one of the p-th roots of A: = A, moreover, in this case the integral does not depend on 
Tfc with k divisible by p |10j. For non-monomial potentials V{x) there are essentially different choices 
of L and essentially different Kontsevich phases (in their simplest phase, non-monomial potentials are 
reduced to the monomial ones, see [T8| [50]). 

• As usual for matrix models, the original integral (jl.2p is not an adequate definition of the partition 
function: as it is, it describes reasonably only some of the phases. At the next step, it should be 
substituted with a set of differential equations w.r.t. the time variables t^ and Sk so that the partition 
function is defined to be a generic multi-branch solution to this system, with (|1.2p providing integral 
representations for some of the branches. These equations have a simple form of continuous Virasoro 
constraints for the simplest case of V{x) = [SSllMllMlIll], i.e. for the original Kontsevich model [9], 
become more sophisticated M^^^+^^-constraints for the GKM with monomial potential [n\ [2^ [571 13S] 
and turn into even more sophisticated relations for a generic V{x), especially when s-dependence 
is also taken into account. For fixed V{x) the different branches in Kontsevich phase possess loop 
expansions and are further associated with shifts t^ — > Tj. + t]^^ so that expansions are in positive 
powers of t- variables with T appearing in denominators - just like in the case of the Hermitian matrix 
model. The (p, g)-model is the GKM with Vp{x) = ^q^j- in the phase with T/; 7^ for k = \,p + q only 
[23J. 



• One of the most remarkable properties of GKM is the p-q duality [T71 [181 II]- the relation 
between partition functions Zp^q and Zq^p. It is not a literal coincidence between the two branches of 
the partition function, instead they are associated with two different coverings of one and the same 
spectral curve and should coincide after an appropriate change of time- variables. 



^Note that p here can be negative equally well, the anti-polynomial Kontsevich case, |34] 
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• Since partition functions of all models, associated with the Itzykson-Zuber integral (jl.2p possess 
the determinant representations ()1.4p . it is natural that they are r-functions of the KP and Toda 
families in t- variables [TOl [HI E] • They also possess certain integrability properties w.r.t. the s- 
variables [181 El EOl [50]. Thus, the entire theory of Itzykson-Zuber-Kontsevich models is indeed a 
piece of theory of stringy r-functions. These r-functions are, in fact, closely related to the Hermitian 
model r-functions: both classes belong to the same matrix model M-theory [5]. 

• The main questions to be addressed in the course of study of every particular phase of every 
particular model are listed in the following table: 



IZK integral = 
partition function / = exp g'^^~'^ F^^^ 

i 

Ward identities in the form 
of Virasoro and W constraints L{z)I = 



^enus-zero part F^'^'^ of partition function 
i 

bare spectral curve Eg 
i 

full spectral curve S, 
describing the matrix-model r-function 



multi-resolvents p''P^^\zi, . . . , Zm) 
i 

poly-differentials on Eg [ 

i 

Laplace transforms 7/(pI™)(xi, . . . ,Xr, 



The bare spectral curve Sq is an important characteristic of the branch of the partition function: 
different phases of the same model differ by the shape of Sq. In fact, in order to describe higher 
multi-resolvents in generic phases [Ullj, the spectral curve should be made dependent on the genus 
expansion parameter (7, thus breaking the simple association between Sq and F^^'^ . The study of this 
phenomenon can be one of the clues to the construction of the last vertical arrow in the left column, 
relating Sq with the full spectral curve S. As every KP/Toda r-function, the partition function is 
formally associated with a point of the Grassmannian [39j and, thus, formally with some infinite 
genus Riemann surface: this is exactly what we call S. The horizontal line in the center of the 
table is a functorial map from complex curves to a hierarchical family of poly-differentials, which can 
be described and studied independently of other parts of the table. A big step in this direction is 
described in [12], but representation in terms of free fields on Sq is still lacking, even in the simplest 
phases. Moreover, this map depends on additional conditions imposed on the poly-differentials, which 
are actually related to the choice between different branches of partition function, made also beyond 
the genus zero approximation. Also lacking is a description of the vertical arrow from poly-differentials 
to their Laplace transforms, which should be very interesting, because the Laplace transforms possess 
a very simple m-fold integral representation, at least, in some phases jl6j . 



2 Kontsevich model 

2.1 Solving loop equations 

2.1.1 Kontsevich model: definitions 

We start with defining the Kontsevich model. As explained in the Introduction, we define any matrix 
model partition function as a solution to an infinite set of equations. In particular, the Kontsevich 
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partition function is defined to satisfy the continuous Virasoro constraints: 

L^{z)Z{t) = 0, 

Ln{dzf 



L.{z) = [:.P{z)-)_=Y, 



z 

n 



n+2 



1 °° r 

fc=0 ^ 



d 



(2.1) 

In order to define the branch of the partition function, we shift the times, 

t2k+i = T2k+i + T2k+i 0<k<N, (2.2) 
and consider the partition function to be a formal power series in the shifted times T2k+i that satisfies 

Now we shah fohow the hne of paper [Ij and rewrite ()2.ip in the form of loop equations that admit 
recursion solving. To this end, we introduce 

• the loop operator 

^/ ^ \ ^ 1 d 



z' 

n>0 



+3/2 dT2n+l (2-^) 



the generating function for Tf^ (a polynomial of degree N) 

N 

W{z) = Y.i'^k + l)z''~^I^T2k+^ (2.4) 

A;=0 

the generating function for (a power series) 



v{z) = J^(2A; + l)z^-i/V2fc+i (2.5) 

fc=0 

• the projector onto the negative part of series 

P- I z'^aX = ^'«A. Pt = ^- Pz (2.6) 

lfc=-oo J k<-l 

• the free energy and its topological expansion w.r.t. to the genus p 

Z{t) = e^^^^) T = Y 9''^^^ -^1^=0 = F[T] (2.7) 

• the generating resolvent 

G{z\t)=V{z)T{t) (2.8) 

• and the multi-resolvents 

pl™)(zi,...,z™) =V(zi)---V(z^).F|,=o (2.9) 
p(Pl-)(zi, . . . , z„) = V(zi) • • • V{z^)T'^p\=o (2.10) 



• the set of /-functions generated by the i?-check operator 

P- {W{z)G{z\t)) = W{z)G{z\t) - /(z|r) (2.11) 

f{z\T) = P+{W{z)G{z\t)) = R{z)m (2.12) 

N-2 N r. 

^(^) = E E ^'~'"-'(2fc + l)W— (2.13) 

/(^'l"^+i)(z|zi, . . . , z„) = V(zi) • • • V{zMz)^^''\r)\r=o = R{z)p'^P\^\z^, ...,zm) (2-14) 



2.1.2 The loop equation and recursion relations on the multi-resolvents 

Now rewrite the Virasoro constraints ()2.ip in the form of the loop equation for the resolvent 



P;{v{z)G{z)) + W{z)G{z) - f{z) + \g\z) + ^\{z)G{z) + ^ + ^^^^ = ° ^^'^^^ 
Applying the operator V to this equation k times, using the identity 

V{X)P- {v{z)h{z)} = 2d^ { ^^^^ } + {v{z)V{x)h{z)} (2.16) 

and ultimately putting all Tf. = 0, one comes to the set of recursion relations for the multi-resolvents, 



i=l 



2 

k\+k2=k 



-.2 

■ :nt,.r(2-fc)!^ 

These recursive relations are invariant with respect to two different scaling transformations with the 
following scaling exponents: 



degT2n+i = n-l deg5 =degF = -3 degr2„+i = n-l 

5 5 (2.18) 

degz = -l degV = - degp^^l^) = -3 + -fc 3p 

and 

deg' g = deg' U = I deg' F = 2 

(2.19) 

deg' z = deg' V = -1 deg' p^^'l''^ = 2 - A: - 2p 
Making the genus expansion of the recursive equations (j2.17p , one obtains for the g^^-iexin 

2y:4. { ''>'^""'<"'--^--^-'-^""'^' -''>j-/H^^.)(.|......,.,H 

Z Zi 



+H^(z)p(p|'^+i) (^, ^,) + i ^ /'^'"^^'^ (^, , • • • , )p(^-^i'=^+^) (^, . . . , )+ 



g=0 ki+k2=k 



^ 1 ^(p-i|A.+2) , . . . ^ + l-5k,oSp,i + ^ 3/2(2-fc)i ^i"''^P-° = (2.20) 

These double recursion relations (in p and k) can be used to determine all the multi-resolvents 
recursively. E.g., for we have a quadratic equation 

(p(°|i)(z))2 + 2W{z) ■ p(°|i)(z) - 2i?(z)F(0)[T] + rf/z (2.21) 

its solution bein^ 

p^^\'^\z) = -W{z)+y{z) (2.22) 
where y{z) is a multi- valued function of z 

y2 = {W{z) f - ^ + 2i?(z)FW [T] (2.23) 

Making further iterations, the multi-density ^(^1^^+^) enters (|2.20p linearly with the factor y{z), and 
one can make iterations, e.g., using some computer algebra system in the following order: (0|1) — > 
(1|0) - (0|2) - (1|1) ^ (2|0) - (0|3) - (1|2) - • • • 

Note that the recursion relations contain a lot of ambiguity encoded in the functions f^P^^\z) = 
R{z)F^P')[T]. Indeed, F[T] can be an arbitrary function that satisfies the two constraints L-i, Lq 
(these are and Lg-constraints with all Tfc = 0): 

g(2*+l)r..«_ + ^ = gp.+ l)T.„«^ + |=0 (2.24) 

Therefore, the space of solutions to the loop equations (Virasoro constraints) is parameterized by such 
functions F[r]. 

2.2 Solving the reduced Virasoro constraints 

The general solution of the second equation of (12.24p . i.e. of the Lo-constraint, is 

F[T] = - 8(2^+1) ^°g^2^+i + ^(^1' • • • ' (2.25) 

where 

_ (2Ar-2fc + l)!! -2T2N-2k+i 

to AT I IMI " 2]V-2fc+l (2.26) 

(2iV + l)!! ^_2j.^^_^^^^^ 



and F is an arbitrary function. Then the L_i-constraint reads as 

dxk+i 



Vv - m+m\ ^.2 .2 271 



k=0 

where xo = 1 is not an independent variable. Its general solution is 
-:, [(2iV + l)!!]2 r^i 



XN(.m,m m)dfji + F{r]2, . . . , w) (2.28) 







where F is a new arbitrary function, and we made the triangle change of variable^ generated by the 
following relations between the generating functions 

N N 00 , sp 

X{z) ^ ^(^) = E ^'^^^ = E ^X^(^) (2-29) 

k=l k=l p=l ^ 



^Note that the non-meromorphic term W{z)dz in p^-^^^^dz just cancels the singular part of the z-expansion of ydz at 
infinity. 

^Since this change of variables is triangle, the variables x* can be equally well expressed through ry, . 
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In other words, r]k for k < N -\- 1 are defined from the equation e'^^'^^ = ~^ 0{z^'^^) with all Xk 

and % equal to zero if A; > A^. 

In order to prove that (j2.28p . indeed, solves (j2.27p . one suffices to note that, since 



the L_i-constraint (|2.27|) reads as 

^1 OF ^1 % dF ^^Vr W dP £ dx{z) x^-\z) ^ 

dF^ y f x''^\z)[l + x{z)] ^_ ^dF f]_^^^_dF_ 
i^'li^i Jo ^ ^dr]i%zi ^ dm 



(2.30) 



Note that this choice of r/- variables is in no way unique: e.g., one can equally well transform 
r]2, - ■ ■ ,'r]N to any new N — 2 variables without changing the formulas above. For instance, one can 
request that the transformation is linear in all xi with / > 3, 

k 

m = Qk{xi,X2) + xiPkAxi,X2) (2.31) 

l>3 

and check that there exist polynomials Qk and Pk^i that preserve relation (I2.30p . Inserting ()2.3ip into 
(j2.30p . one immediately obtains that these polynomials satisfy only the equations 

-3^ = -X2-P/c,3 ~H7^ = -Pk,i+i (2.32) 
where we changed xi,2 for new variables ^1^2 

^2 

6 = Xi 6 = X2 - Y (2-33) 

Equations (j2.32p also have a lot of solutions. In particular, one can add to an arbitrary function 
of ^2- Choosing, e.g., this function to be zero and Pk,k = 1) one immediately obtains 



tfc-2 ck 

A; -2^ 2fc(A;-3)! 



(2.34) 



Pk,k-i = {-y^ Qk = {-] 
2.3 Example of transition: = 2 ^ = 1 

Let us consider the limit of ^ and see how the space of solutions to the loops equations, which 
is parameterized by a function of one variable in the case of = 2, reduces to the only solution in 
the case of A^ = 1. 
For N = 2 one has 

- 1 -2^3 _ 1 -2Ti 

" 5 (-2T5)3/5 ~ 15 (-2r5)V5 (2-^^^ 

Vi = -Xi m = ]pd - X2 (2-36) 

F[T] = -|! logTs - 152 (^x:i - ^XIX2 + \xixi) + hm) (2.37) 
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In order to have a smooth transition as T5 goes to zero, one has to cancel the singularity, i.e. F{r]2) 
must have an asymptotics 

2 

F(r/2) = 15^/2 -T?^/^ -f- log 772 + 0(1), m^oo (2.38) 
and, therefore, one obtains 

F[r](^=i) = hm F[r](^=2) = -^logri--^^ + const (2.39) 

which coincides with (j3.2p below. 

One may use the scaling symmetry ()2.18p in order to further restrict the function F. Indeed, using 
that the scaling dimension of r/2 is —6/5, one immediately comes to the expansion 

00 

^ = E /"Cp^f '^'""^ + C[9^ log m (2.40) 

if assuming that the symmetry (|2.18p does not change the solution, or, putting this differently, that 

F does not contain any additional dimensional parameters changing under this symmetry transforma- 
tion. Co and C( can be obtained by comparing (j2.40p with (j2.38p . Therefore, the asymptotic (j2.38p 
corresponds actually to the semiclassical limit of — > 0. 

Geometrically this transition corresponds to degeneration of the torus into the sphere. Indeed, 
formula (|2.23p describes the torus in the case of = 2. As it follows from (j2.37p and (|2.38p . 
2i?(z)F(°)[T] = lQT^dF^^\T]/dTi ^ as T5 ^ 0. Therefore, ((2:23]) transforms under this tran- 
sition exactly into the sphere corresponding to = 1. 



2.4 Resolvents 

As as explained above, starting from the one-point resolvent 

p(0|i)(z) = -w{z) + y{z) 
one recursively calculates further resolvents: 



(2.41) 



^{0|2) 



(3^1,2:2) 



y{xi) 



2d, 



X2 



p(0|l)(xi)-p(0|l)(x2) 



Xi - X2 



/(°I^)(X1|X2) + 



3/2 



(x) = ^f>|2)(x,rr)-/Wi)( 



QX 



(2.42) 



(2.43) 



etc. 



These resolvents can be given a geometric meaning. Indeed, the r.h.s. of (j2.23p is actually a 
polynomial of degree 2A^ — 1. Thus, equation ()2.23p defines a hyperelliptic curve Sq = C of genus 

— 1 in a generic case. (This genus should not be confused with the genus corresponding to the 
expansion of the free energy in powers of usually labeled by p.) This bare spectral curve is actually 
essential for constructing all the multi-resolvents because these are meromorphic multi-differentials 



P 



(p|m) 



P 



(p\m) 



i)dzi---dzm on this curve with specified singularities {p^^^^^ and p(^l^) are 



distinguished differentials playing a specific role). Typically this leaves some room for adding holo- 
morphic differentials which exactly corresponds to the ambiguity in solutions to the loop equations. 

Indeed, because of equations ()2.24p . one has an arbitrary function F of {N + 1) — 2 = A^ — 1 
variables, and, at each step of recursive computation of the multi-resolvents, one finds in p^^l™-) some 



ambiguous terms 



Qm p{p) 



< TLi < N — 2. Fixing these terms is equivalent to fixing the 
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periods of All these terms are certainly fully determined by 



dF(p'> 

dT2n + l 



i.e. by fixing the periods 



of 

As an example, consider the two-point resolvent. It can be rewritten in the form 



,(0|2) _ (0|2) (0|2) _ (0|2) 
^ ~ rhol ^ rglob yioc 



where p~^^^ is a holomorphic bi-differential on the curve C, 

(012), s k{xi,X2)F'^^\T]dxidx2 
y{xi)y{x2) 



(2.44) 



(2.45) 



Pgiob ^ meromorphic bi-differential on the curve C that has the singularity at xi = X2 at the both 
sheets of the curve of the following type: 



(o|2) 2dxidx2 , 



(2.46) 



At last, pL?^ 



„(0|2) 
Pglob 



has the same behaviour as p~„i\,j^ at infinity and cancels the singularity 

y{z)=y/z 

(I2.46P of p'gi^il at one of the sheets of the curve. 
Manifestly, 

(012) [xi + X2] dxidx2 

' (2.47) 



J0|2) 



Pi. 



Pflobi^i^^^) 



where 



y{xi)y{x2){xi - X2] 

B{xi,X2) ■- 

C(xi,X2) = 
N-2N'-2 N N 



^2'^ (^1 ~ ^2)^ 

2 ((Xi + X2)S(X1, X2) + (Xi - X2fC{xx,X2) + i(xi, X2)FW [T]) 

(2.48) 



W(xx)W(x2) 



+ 



W{X2) ^ W{xx) 



X1^X2_3 
X2 X\ 



X\X2 



k—ra—2k'—m'—2 



x\ 



i2k + l)i2k' + l)T2k+iT2k'+i 



d 



Kix„x2) = Y.E E E -1 

m=0 m'=0 k=ni+2 k'=m'+2 
N k-4 

A{xi,X2) = - X^Z] {^^""^[3^2 - 5x1 + 2n(x2 - xi)] + (x2 ^ xi)} {2k + l)T2k+i 



d 



dT2m+l dT2m'+l 

d 



k=4: n=0 



+ 



2n+l 



+4 



d 8 7 d 

5T5- — h 7T7- h -(Xi X2)r7^— 

^dTi ^ dTi 2^ ' ^ dTi 



Note that the numerator of p^gi^^ is actually a polynomial in xi, X2. 

Since is a holomorphic bi-differential, the second derivatives gy^^^^y] .^^ (entering this dif- 

ferential) control the periods of ^('^l^) and do not affect its singularities. Note that the number of 
independent variables in F\T] is equal to the genus of C: — 1. 

Similarly one can deal with other resolvents in order to check that the multi-resolvents = 
p(p|™)(2^^ • • • , Zm)dzi ■ ■ ■ dzm are meromorphic multi-differentials (except for the cases (0|1) and (0|2)) 
on the curve C and generically have poles of order 6p-|-2m — 4 in points (and only in these point) where 
y = 0. These multi-resolvents can be further restricted with using symmetries (j2.18p and (j2.19p . 

We discuss these general properties of multi-resolvents and its applications in more details in the 
next section in the simplest example of the Gaussian Kontsevich model. 
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2.5 CFT representation 

As already mentioned in 1 
agator in a certain CFT: 



As already mentioned in the Introduction, the two-point function p^g}^^ can be represented as a prop- 



P^Xh^i^^2) = {dX{z,)dX{z2)) (2.49) 

where X is some local field defined on CP^ parameterized by z. This is because there is a singularity at 
zi = Z2, i.e. when the arguments of the fields coincide. The CFT is defined by the covering C CP^. 
One way is to consider the scalar field living on the C as a collection of two fields Xi , X2 living on the 

corresponding sheets of the covering: C C . Then there will be a monodromy Xi ^ X2 when 

z goes around a branch point y = 0. Then the field X is a linear combination of fields Xi, X2 that 
diagonalizes this monodromy: X = Xi — X2, X ^ —X. 

To put this differently, let z parameterize the whole world-sheet now (i.e. topologically it would 

be a sphere), but the target space is now an orbifold: C/Z2 — CP""^ C/Z2. The branching points 
y = are now just points where string wraps around the Z2-fixed point. 

Both of these approaches can be actually described in the same way via the branching point 
operators of 001 HI]- To this end, one needs to introduce the twist field ai/2iw,w) with the following 
operator product expansion (OPE): 

dX{z)ay2{w,w) ~ (z - w)'^^\y 2(^,^0) + ... (2.50) 

where T1/2 is sometimes called excited twist field. Then, one can write 

{dX{zi)dX{z2)) = /dXizi)dXiz2) H ai/2{m,Wi)\ (2.51) 

\ y(wi) = I Q 

By (•)q we denote the correlator in the ordinary CFT on CP^. This would provide us with the necessary 
structure of singularities 

{dX{zi)dX{z2)) ^ ^—2+0(1), zi^Z2 (2.52) 

[zi - Z2Y 

{dX{zi)dX{z2)) ~ {z, - Wj)-^l^, z, ^ W.J, y{wj) = (2.53) 

from which one can deduce (j2.49p . 

One can also include in this correlator. It would control the global monodromy properties of 
the filed X, i.e. how it changes when z goes around the cycles on C. 

At last, 

pf}^\zuz2) = {dX{z,)dX{z2)ay2{m), (2.54) 
since p^l^^^\zi^ Z2) knows nothing about the branching points Wj, see (|2.47p . 

3 Gaussian branch of Kontsevich model 
3.1 Specific of the Gaussian branch 

In this section we consider the special, simplest case with only the first two times non-perturbatively 
turned on {N = 1). Then, the i2-check operator ()2.13p identically vanishes 

%) = ^ /(^'l'^) = 0, yp,k (3.1) 

Therefore, there are no ambiguities in resolvents in this case. This key feature suggests a separate study 
of this distinguished case. This case is also a counterpart of the Gaussian branch of the Hermitian 
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matrix model, hence, we call it the Gaussian branch. To avoid misunderstanding, note that it has 
nothing to do with the Gaussian integral! 

Given Ti = a and T3 = — 5M, the solution of (|2.24p is 

F[T]= F{M,a) = -^- ^ log ^ (3.2) 



and the curve is 



y' = MHz - ») (s = 1^ j (34) 

Therefore, the resolvents non-trivially depend only on one parameter s (M can be effectively removed 
by rescalings). To simplify formulas, we consider from now on the redefined curve 



Y{z) = y{z)/M = ^/Y^s (3.6) 

In the Gaussian case, the recurrent relation (j2.20p can be simplified. More concretely, for suffi- 
ciently large indices it can be written in the form 

y (.)p(^l'=+^) (., .1 , . . . , z,) = -2 ^ a, P!^!H£klllli^+ 

VI +P2 =P 
fci + Pi > 

(3.7) 

where the subscript (•)mcr means that one has to replace p^^^"^") with p'gi^^ leaving all other p's unchanged. 

Indeed, for large enough indices ()2.20p contains, in the Gaussian case, only four terms. The only 
term containing pC^I^) combines with that containing W{z) to produce y{z), while p'l^}^^ in the sum 

quadratic in /o's cancels the non-meromorphic (~ (^)^) part of the first term. 

Recurrent relations ()3.7p celebrate an important property that leads to drastic simplifications in 
the Gaussian case, which allows one to get rid of the only parameter s: 

Important formula: All the s- dependence of the resolvents (except for p = 0, m = 1 and p = 0, 

m = 2 cases) is actually encoded only in the differences Zi — s: 



P 



(P|™)(Z1, . . . , Zn,\s) = />(^'I™)(Z1 -s,...,Zm- s\0) (3.^ 



In order to prove this formula, let us denote though Li the first-order part of the differential 

operator /_i = Yl (2^ + l)i2fc+i gj^-^) -^-i = ^-i + Then, the L_i-constraint on T is 

k=l 

1-i:F=-^ (3.9) 

To prove (13. Sp . one suffices to note that 

[/_i,V(z)] =2— V(z) and /.^l^^^ = sTg— = -2— (3.10) 
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and, using these formulae, to show immediately that — ^ and ^ ^ acting on 

i 

. . . , z^) =^ V(zi) • • • V(z„)^(f)|.=o (3.11) 

are equal to each other whenever (p, m) ^ (0, 1) or (0, 2) (so that one can ignore the r.h.s. of (j3.9p ). 

Formula ()3.8p can be also proved by induction using the recursive relations (|3.7p . Indeed, the 
claim is correct for p^^^ by an immediate check. Further, if all the p's in the r.h.s. of (j3.7p enjoy 
the property (|3.8p (by the induction assumption), this is also true for z — zi and y[z) and, thus, for 
p{v\k+i) in the l.h.s. 

We can use y as a standard coordinate on our CP"*^, since z = 

y2 + s. Then says that the 

densities can be written in terms of Y^s only. Thus, the case of arbitrary s is, in a sense, equivalent 
to the s = case. 

3.2 Resolvents 

In the next subsection, we present for a reference manifest expressions for several first densities. They 
all can be obtained recursively using (j3.7p . by hands or with the help of computer (using, e.g., maple). 

3.2.1 First resolvents 

The one-point resolvents: 

• Genus p = 



p(°|i)(z|s) = M 

• Genus p = 1 

• Genus p = 2 

• Genus p = 3 

• Genus p = A 



r- s n\ ..s^ s'^ r(n + l/2) ^ ^ 

Vi-'^^'V =^2g .n.VMn;i)!r(l/2) (^-12) 



.oHN, , , 25025 1 1 , ^ 



.4ll^, , , 56581525 1 1 , ^ 



The list of resolvents grouped by the genus p (we put s = in all resolvents but p^''!^) and /O^'^l^^; 
the s-dependence can be easily restored using formula (|3.8p ): 

• Genus p = 

p^^\^\z\s) = M [-^z + ^ - Y{z)^ (3.17) 

(012)/ I X Zi+Z2-2s ZI+Z2 



(3.19) 
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.3|.) = ^ ^3(,^)^3|,^)y3(,3) (3-20) 



3 1 

p(°l^)(2;i, Z2, Z-i, Zj^\s = 0) = ^2 5/2 5/2 5/2 5/2 ' ^2 ^3 + ^^l 2=2 + ^^l 2=3 2^4 + 2^2 ^3 Z4) 



^1 -^2 -^3 -^4 



(3.21) 



p(°l^)(zi, 2:2, 2^3, 2^4, ^5|0) = j^^3 ^^7/2^^7/2^37/2^^7/2^^7/2 (6^l'^5^4'23'^2 + 5 Zi222'z5223' 

+ 5 zi' Z2' z-i' z/^ + 6 2:1^-2:2^^5^-24-^3 + 6 zx' Z2Z^ z^ zj, + 6 .21^.25^-24-2:3^-22 + 5 zi'z^z^Z'^ 

+6 Zi'^ Z2^ Za,'^ Z4Z5 + 6 Zi'^ Z2^ Z5Z4^ Z3 + b Zi' Z2^ Z^"^ Zi' + Q ZiZ2^ Zr,Z4^ Zz' + () ZiZ2Z5^ Z4^ Z3^ 



+6 ZiZ2'^ Z^'^ Z4Z2,'^ -\- ^ Z\Z^ Z/i^ Z'iZ2' + h Z2^ Z^'^ Z^^ Z/^') 



(3.22) 



Genus p = 1 



(3.24) 

^1 ^2 

?W3)(zi, Z2, 23IO) = • (35 Zi^ 22^ + 30 Zi^ ^2^ Z3 + 30^1^^2232 + 35 Zi^ 2:33 

-2l •22 2:3 

+30 Zi^ Z2^ 23 + 18 Zi^ Z2^ ^3^ + 30 Zi^ Z2 ^3^ + 30 2i Z2^ Z^^ + 30 2i 22^ ^3^ + 35 Z2^ Z^^) 

(3.25) 



p(l|^)(2:i,Z2,Z3,24|0) = 8 ^4 ,^11/2,^11/2,311/2^^11/2 (105 21^2^324^ 
+ 100zi^Z2^23^24^ + I05 2i^z/Z3^ + 90^1322^23^^4^ + 105 23^22^2421^ + 5423^223243213 
+ 10021^22^23^24^ + 9023^22324321'' + 90 23322324^21'' + 100 2i'' 22^23^24^ + Z2^ Z^^ Z^^ Zi^ 
+ 105 23 ''22^24213 + 9022^3324^213 + 105 23322'' 2421^ + 9023^23324321^ + 90 23^22324^213 
+9021322223324"^ + 10521^24^23223 + 90 21^22^233243 + 90 22^^233243212 + 105 22''23324''2i 
+9022^232243213 + 10022^23224^212 + 10022^23^242212 + 10022^23224221^ + 105 21324^23^22 
+ 105 21^24323'' 22 + 10522^23^24321 + 105 22 Si ''23" + 105 22323 ^^24^1 + 90 22323324''2i2 
+ 105 22^23^24^^ + 105 21^^24^23322 + 105 22^21'' 24^^ + 105 22^2324^^213) 

(3.26) 

Genus p = 2 

.011 w , , 105 1 1 , , 

p(2|i)(^|g) = ,,3„iw . (3.27) 



128 M3 yii(2) 
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or 1 

P^'^'Hzi, Z2\0) = ,,3^4 13/2 13/2 • (33-1^ + 27.1^.2 

^1 ^2 

+29 zi^ Z2^ + 29 zi^ Z2^ + 27 zi + 33 Z2^) 
35 

p(^|3)(z„Z„Z3|0) = ^5 ,^15/2,^15/2,315/2 (396z,^Z,z/ 

+396 Z2^23^^i + 396 zx^Z2^Z'i'' + 406 Z2^zx^zi' + 396 Z2^^i^23^ + 396 Z2^z^Z'i 
+396 ^2^21*^23 + 348 zx^Z2^zi + 324 z^^ Z2^ z.^ + 396 z^^ Z2^ z^ + 324 zs^za^zi^ 
+348 z-i'^z-^z^ + 396 za^za^'-Zi^ + 429 Z2^zz^ + 348 23^22^21^ + 396 zi^azs^ (3-28) 
+396 Z2^zx^Z'i^ + 324 ^2^21^23^ + 406 ^1^22^23^ + 348 zx'Z2^z-i^ + 429 21^3^ 
+348 21^23^2^ + 396 22^21^23^ + 360 ^2^21^3^ + 406 zx^Z2^z^^ 
+ 396 Z2^zx^zz^ + 348 ^2^21^23^ + 429 zx^Z2^) 

• Genus p = 3 

, , 25025 1 1 , ^ 



35 

p(3|2)(,^,,2|0) = 1024 ,^19/2,^19/2 (^^^^^^1' + 10725 Zi^Z2 + 11011 Za'^l' + 11066 21^22 



+10926 zi'' 22^ + 11066 22^21^ + 1101121^22^ + 10725 22 '^zi + 12155^2^) 



3.2.2 Resolvents: general formulae and relations 



(3.30) 



For generic p and m, (p, m) 7^ (0, 1) or (0, 2), in the Gaussian case all the resolvents are of the following 
form: 



(2;i,...,2:m|0) - ^2p+m-2 J-J™ ^ y6p+2m-3(,.) J (3-31) 

where Qp^m are homogeneous symmetric polynomials in {zj} of degree degQp_m = {ra— l)(m + 3p — 3). 
For m = 1 degQp^m = 0, and Qp.m is just a constant 

so that the formal power series p'^H-^) — X^^o ^^''/'^^'^H-^) can be converted into the hyper geometric 
function. Its z-expansion is 



(6p-3)!! 1 1 ^g" r(n + 3p-l/2) 



23P3P„! M2P-1 ^3p-l/2 ,n ^1 r(3i5 - 1/2) 

which gives the general formula for the ((o"o)"(T3p_2+n) intersection numbers [9]. 
There is also a general formula for the genus zero [p = 0) multi-resolvents: 
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A general formula for the p = 1 multi-resolvents is less explicit. It can be written with the help of the 
generating resolvent ()2.8p treated as a functional of (p{z) such that d(p{z)/dz = (j2.5p : 



A:=0 i=l'^^ lli=l-^ ^^J 



(3.35) 



fc=0 



C - s)3/2 



s=0 



<j>(z)dz 
'C (^2_s)3/2 



s=0 



z" — 2t) and the contour C encircles oo. Then, the logarithm of the generating 



where = (z^ 
resolvent generates the p = 1 multi-resolvents 



fc=0 ■ i=l 



(3.36) 



i.e. the genus one multi-resolvents are connected parts (up to the factor of 24) of the /c-point function 
generated by the generating resolvent, 



p(^l^)(-i,. 



1 (51ogGW{0} 



24 



2:i)...5(/>(zfc) 



(/)=0 



For example. 



pWi)(.) 



1 (5G(o){0} 



1 



24M2 



^(l|2), 



24 



2:1,^2) 



1 



24M 



fj,{z)dt 



24M 



Kz)\t=o 



8Mz5/2 



1 

24 



,52^(0) {0} _ SG^^^c/)} (^GW{(/>} \ 



zi)(50(2;2) 5(^(zi) 



0=0 



24M2 I 7 



/ fi{zi)fi{z2)dt f ^{zi)dt f fi{z2)dt 



2[fl{zi)fJ,{z2) + 2fl{zi)il{z2) + ^i.{Zl)|l{z2)) 



t2 J t2 



fl{zi)fl{z2) 



5zf + 3ziZ2 + 5^2 



t=0 



8M2z; 



7/2 7/2 



3.2.3 Proofs and comments 

In the forthcoming considerations, we use the symmetries (|2.18l ^ and ()2.19p which should be sup- 
plemented with weights of genus-expanded multi-resolvents under the first transformation, (j2.18p : 
deg p^pI*"') = —3 + + 3p, and under the second transformation, ()2.19p : deg' p^P^^^ = 2 — k — 2p. 

All the proofs in this subsection are done by induction. Note that we often omit as trivial checking 
the induction base. An equivalent way to obtain multi-resolvents is described in s l4.3l 



*For Al — 1, s = this symmetry requirements are equivalent to the standard claim that the free energy has an 
expansion in r such that 

fe 



rii 



3p-3 



(3.37) 



C's being some numerical constants. This is usually derived [9] from the fact that the intersection number of a collection 
of forms is non-zero only when the sum of their degrees is equal to the dimension of the manifold. 
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Proof of (j3.3ip : Using (j3.7p it is easy to show by induction that p^^^'^\{zi}) is meromorphic on 
the CP^ in each Zi and can have poles only at y = 0. Then, in order to prove (jS.Sip . it is enough to 
show (due to the symmetricity in {zi}) that for some Zj the order of pole in 1" = is (6p + 2m — 3). 
This also can be done straightforwardly by induction. Then the degree of the polynomial Qp^m is 
completely fixed by symmetry (|2.18|) : deg Qp^m = (m — l)(m + 3p — 3), while the power of M in (j3.3ip 
can be determined by symmetry ()2.19p . □ 

In fact, the second symmetry (|2.19p allows us to put hereafter M = 1, while the first symmetry 
(I2T8I1 does not affect M at all. 

Proof of (j3.34p : Without any loss of generality one can put s = (i.e. y{z) = z^/^) using formula 
(j3.8p and identify ds = — Yli=o ^Zi- From the definition of p's and using symmetry (j2.18p . one obtains 



P H{^^}.=l) - 2^ ^ (3.38) 

In each term of the sum at the r.h.s. of this expression there exists i such that rii = 0. Let us denote 
through the subscript A the coefficient in front of in the asymptotic expansion at z ^ oo. Since 
p is symmetric in Zj's, 

Pa {{^^\^=l)= 2^ U-1 (3-39) 

contains the same information as p^^^^^ itself. Thus, p can be in principle restored from pA- However, 
we do not need to do this explicitly. Just assume ()3.34p is correct for k < K (\t is trivially correct for 
K = 2). To prove it for k = {K + 1), it is enough to show that 

pT^'\{z.}1.) = f-2Ea.') (3.40) 



nr=i 



One can easily do it by considering ~ ^ term of the asymptotic of the recursive relation (13. 7p : 

pT^'\{z^}f=,) = [-^jZd^ P^°I^H{-af=i) (3.41) 
Then, (j3.34p is correct by induction. □ 

Proof of (j3.36p : Similarly to the genus zero case, put s = and write 

p^ >{{zi}i^^) - ^ — (3 42) 

Y.ni=k lli=l 

Now either there exists such i that = or Vi nj = 1. Therefore, in contrast with the genus zero 
case, all the information about 

pirn is 

contained both in 

Pa {{^^}^=l)= 2^ 3/2+n. (3.43) 

Y,ni=k llj=l 

and in Ck = 3^'^Ci...i . First, we deal with this c^. One can easily construct such that §(j dx^^^j^ = 
^(5„.o- Now, we are interested only in the terms where each dg acts once on each ^^^-p- in the k-point 
function in (j3.35p . Then, (j3.35p reads as 

GM=E^-^' = rT^ (3.44) 

fc=0 
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and 

oo ^ 

log G[4>\ = - log(l - a) = ^ ^ (3.45) 



k 

k=l 



i.e., (|3.36p is equivalent to the equality 

(k-iy. 



Ck = 



24 



(3.46) 



Now, as usual, we prove it by induction. Assume that, indeed, Ck = {k— l)!/24 for k < K. Computing 
the term — of the asymptotic of (j3.7p . one obtains (again only the first term at the r.h.s. 



survives) 



= -2K{3^ck) + 5K{3''ck) = 3^+^Kck (3-47) 
where the two terms come from differentiating the denominator {z — Zi) and the numerator respectively 
(one then expands ^^5^. = \dzi + ^Zidz^ + ■■■)■ Therefore, ck = ^^24^^' ■ 
As the second step, we study p^'^''- The ^-asymptotic of the (|3.7p gives us 

pT^'\{^^]U) = (-2E^^) P^'^'\{-^)U. (3.48) 

Introduce the notation p'^^\{zi}\^^) = {2ds)^Y^ — ^ 3^ ^'^'^ p'^^ln{{^i}\=i) for its connected part. Now 

we again apply the induction and assume that, for k < K, p^^^^\-) = -^piolni')- To prove it for 
A; = IT + 1 it is enough to show that 

Due to (j3.48p and to the induction assumption, it is equivalent to 

p;i+:l(.)=2a.pW(.) (3.50) 

For the complete functions such a relation is obvious from the definition: p^^^^\-) = 2dsP^^\-) ■ Let 
I = 1 . . . k and zq = z. The complete {K + l)-point function is expressed through the connected ones 



as follows 



U /,=/u{0}- 

i=i 



E E n p£^.P({-a.e/,)/^ii^'^({^a.e/,u{0}) 

■ / — I T — 1 V * / 



1 = 1 3 = 1 
\J Ij=I J^l 

The i-asymptotic of this equation is 

pT'\{z.ui)= e En/^i^^^(i-'Wp-^;iii:i^({-a../,) 

3 = 1 

Now, using the induction assumption, 

s s 

2a,p(^)({za,e/) = Pcon„,A({^aie/) + E E n Pi^J^({^^}ie/P • 2a,p(ll'l)({z,},,g,, 



(3.52) 



1=1,^. (3-53) 
U ij=i 1 ^ i 
1=1 
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s j=l s>l 1=1 



□ 



3.3 Matrix integral representation 

This partition function Zk can be presented as the Hermitian matrix integral depending on the 
external matrix A, 

[DX exp (-4t^X^ - -^TtAxA 

Zk = ^ (3.56) 

JDX exp (-^TrAX^j 

where the integral is understood as a perturbative power series in T2k+i = g—r- 

— Tr^-2fc-i_ Notg 
2k + 1 

that this integral does not depend on the size of matrices X and A provided it is being considered as 
a function of [lOj. By the shift of the integration variable, it can be also reduced to the form 
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9 N f DX exp(-H^TiX^ + TtAX] 

Zk = exp ( -— TrAi ) ^ ^ (3.57) 

/ J DX expf-^TiAxA 



where 3A = A"^, i.e. T2k+i = Sttt^ Tr A M 1 

3.4 Okounkov's representation of the Laplace transformed resolvents 

Further on in this section we put M = 2, g = 1 and s = 0. 

3.4.1 Laplace transform of the resolvents 

Let us introduce the Laplace-transformed resolvents rj: 

k 



p^''\zi,...,Zk) = 2'' Jfldxie ?''''r/l'=)(xi,...,Xfc) (3.58) 



k 



^ /fJdz,e?"''VI'=)(zi,...,z,,) (3.59) 

i i=i 



(4 

c 

where the contour C encircles 0, beginning and ending at the negative infinity with respect to the 
branch cut along the negative real ray. 

The manifest expressions for a few first r] are 

??(0|3)(xi,X2,X3) = (3.60) 



^In order to introduce an arbitrary shifted first time, t2fc+i = T2k+i wliere M is a parameter, one should 

consider instead of p.56p the integral 

exp(-l|f^TrX='-2g|TrAX^) . ^ , . / OX exp (-ff^TrX^ + TrAx) 
— ^ ^ - oxp TrAs ' ^ ^ 



/L>X exp (-^Tr^A-a) V 6(? J J dX e^p (^-^TrAX 



which is a function of the same T2k+i ~ Q rTrA" ~ — — ^— pTrA" 
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'/ V l'^' o;'^4/ — 2 2^1 2 3 4 ~'*']^ ■*^2 3 4 

, 1/2 3/2 1/2 1/2 , 3/2 1/2 1/2 l/2x 
"I Y "-"2 3 "^4 2 3 4. / 



(3.61) 



r/(^l^)(xO = ^xf (3.62) 

(xi, 2:2) = — ; — {xV"^ x^J'^ + xV'^ + xV"^ xV"^) (3.63) 
24 vr 

, o™l/2 5/2 3/2 3/2 1/2 5/2 5/2 1/2 3/2 3/2 5/2 1/2 (3.64) 
^jli-^ 2 3 ^tij 2 3 / 



77^^l^)(xi, X2) = — {5x\^'^X2^^'^ + 5x\^^'^X2^'^ + 15x'^^^X2^^ 

+15xf ^x^/^ + 29x?/^x^/^ + 29x[/^x^/^) 



(3.66) 



One can easily turn on nonzero s using r/(pl™)(xi, . . . ,Xm\s) = e * r/(pl'")(xi, . . . ,Xm,|0) (except for 
p = 0,m = 1 and p = 0,m = 2 cases). 

The general formula for the genus zero r/-resolvents can be easily obtained from (|3.34|) : 

/ k \ 
E ^i^i 1 / v-^ „ \ 1 



' :^ i=i \ i=i 



lATrj\k J ii 2 I I 3/2 3/2 

^^^^^ I i=l ^ \ i=l ) 2/ ■■■Z^ 

2(47ri)^' Jii'--*- 3/2 _ _ 3/2 2 7r'=/2 ' "^^ ' ' ' 



(3.67) 



3.4.2 One point function 

One can easily verify that the one-point r/-resolvent is as follows 



2 ^ x3/2 

The r.h.s. of p.58p is then equal to 



r?W(x) = -^^ (3.68) 



00 



2 [ dxry(x)e-^^ = ^— [ dxx^P-^/^e-''' 



1 1 r(3p- 1/2) _ ^ 
7^ 2^ nl 19P ^3p-l/2 ~ 2^ 



(3.69) 



a/tt ^ p! 12P 23P-1/2 p! 23p-i 

^ p=0 p=0 



Exactly the same expression is obtained from (j3.33p 

/..f = f;i|_|ll (3.70) 



p!24P22p-i ^3p-i/2 

p=0 p=0 ^ 
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3.4.3 Okounkov's result 



In [T^ Okounkov obtained a representation for the resolvents in terms of finite-dimensional integrals 
which can be thought of as a discrete version of a special functional integral. He introduced the 
functions 



P/ ^ 1 exp(j^^xf) f f ^ {sj - sj+if ^Si + Si+i \ 



and their symmetrized versions 

£^{x) = ^ <S(x<^(i), . . . , (3.72) 

o-eS(n)/{12...n) 

where the summation is over coset representatives modulo the cyclic group generated by the permu- 
tation (12 . . . n). 

Then the resolvent is a generating function for intersection numbers on moduli spaces of curves 
with n fixed points (genus is arbitrary) and is equal (up to some renormalisation and rescaling) to the 
surr[^ 

r?l")(xi,...,x„)= (-l)^(°)+i£:«(x„), (3.73) 
aen„ 

where XI^ is the set of all partitions a of the set {1, . . . ,n} into disjoint union of subsets. For any 
partition a E n„ with £ = £{a) blocks, Xa is the vector of size i formed by sums of Xi over the blocks 
of a. 

For n = 1 formula (|3.73p is very simple: 7]^^\x) = £^{x) = £{x) = and it is exactly what 

we obtained in section [3. 4. 2[ 

Note that Okounkov used the time variables that differ from those typically used in KdV by the 
factor of (2A;-|-1)!! for the (2A;-|-l)th time variabl^. In order to reproduce these factors in the definition 
of the resolvent, one has to use here the Laplace-transformed resolvent r]{x) instead of p{z), since under 
the Laplace transform 

^(2fc + l)!!z'=+3/2 ^ ^ 

3.4.4 L 1 Virasoro constraint and genus-zero resolvents 

In this section we show explicitly that Okounkov's functions (i.e. the r.h.s. of (j3.73p ) satisfy the 
lowest L_i Virasoro constraint. Note that, when proving (j3.34p in sect l3.2.3l we used only the leading 
asymptotics of the recursive relations in z, which is equivalent to the L_i-constraint. In fact, the 
constraint is sufficient to determine the genus zero resolvents in the case when only ^ 0. Hence 
by demonstrating the function satisfies the L_i-constraint one is automatically guaranteed the genus 
zero result is correct (therefore, it is possible just to apply the same arguments as in proof of (I3.34p ). 

The L_i -constraint imposed on the partition function is equivalent to the following constraint on 
the asymptotics of the p-resolvents (compare with (|3.4ip ) : 

pl")(zi, ...,Zn) = -JJ^ pl"-')(^l, . . . , Zn-l) + 0{Z-^/^), Z„ ^ OO (3.75) 

Under the Laplace transform, (j3.59p it leads to 

T7l")(xi, • • • ,Xn) = J— (xi H h Xn-i) T?l""^^(xi, ..,x„_i) + o{y/x^), rc„ ^ (3.76) 



In [16] Okounkov used notation Q instead of 77. 

''We do not care here about some fc-independent factor, because it can be easily eliminated by rescaling. 
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i.e. the r^-resolvents given by formula (j3.73p satisfy the L_i constraint. 
In order to prove this formula, one needs to use the following identity 



oo oo 



da I dbf{a,b) 



_(a-by 

V2Tre 







oo 

dcf{c, c) - ^f{0, 0) + 0(e) ^ y rfc I /(c, c) + ^J^^dJic c)| + 0(e) 

U 

It follows from the decomposition 

/(a, b) = [fia, b) - e(A - (a + 6))/(0, 0)] + /(0, 0) e(A - (a + b)) 



(3.77) 



(3.78) 



where is the Heaviside step function. The difference in the brackets can be proved not to contain 
the 0{y/e) term, while the remaining integral can be computed exactly. 

Now, one can apply to the Laplace transform, (|3.59p of the p-resolvents ()3.75p formula (|3.77p with 

e = 2xk 



£{xi , . . . , Xyi) 



1 exp{^^xl) 



2"7r"/2 n V^i 



'i e 



-Sfc+1- 



i=l 



^."o (3.79) 



exp n S ^ 



2n-l^(n-l)/2 ^ 



JJ dsj ds Ae 



2 



(3.80) 



where 



IXk 



A = l + J^ds + o{^) = 1 

TT 



S - Sk+2 Sk-1 - S ^ Xk-1 + Xk+1 



2xk^ 



2xk-i 



+ o{^) (3.81) 



In S'^ some terms cancels 



in -I) 



Xl^ h Xn-l) 



£'~^{xi, ■ ■ ■ ,X„_l) + o{y/x^) 



(3.82) 



Note that if a is a partition (see [H]), then there are two possibilities: 1) it does not contain the block 
{xn], then 

e'^ixa) = £:^(x„U„=o) + o(^) (3.83) 

2) it contains this block, then 



{l{a)-l)-J^{x^ + ... + Xn-i) 



(3.84) 



In r/ further cancelations take place so that finally one arrives at formula (j3.77p . 

As it was already discussed in the beginning of this subsubsection, one can easily derive using 
(1X751) that 



k \''~^ k 



\i=l / 1=1 

Verification of the Lo-constraint needs much enhanced version of (I3.77P and much more involved 
computations. 
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4 Simplest DV type solution to Kontsevich model: KdV hierarchy 



As we already discussed in sect. 2, there are many solutions to the generic Kontsevich model pa- 
rameterized by an arbitrary function F[T] that satisfies two constraints (j2.24p . in variance with the 
Gaussian case of the previous section, when the solution is unique. However, among all these many 
solutions there is a special family of the so-called Dijkgraaf-Vafa solutions. They are associated with 
a Riemann surface which genus is generically equal to — 1, the number of non-zero times being 
equal to A -|- 1. These solutions are non-generic, and, being extended to depend on higher times T^, 
k > 2 A -|- 1, are required to be associated with the same Riemann surface and to have a smooth limit 
upon bringing these excessive times to zero in order ^2]. Moreover, one typically considers infinitely 
many excessive times (still keeping the genus A — 1 of the curve fixed). Then, F(°)[r] is logarithm 
of the r- function of a Whitham hierarchy w.r.t. these infinitely many times Tj. ]26\ I43| . while the 
complete matrix model partition function as a function of corresponds to a dispersionful integrable 
hierarchy. 

In this section we consider the simplest DV solution, that is, the solution associated with a sphere. 
The partition function Z^it) (as well as Zk[T] = ZK{t)\T=o) of this system is nothing but a r-function 
of the KdV hierarchy, while its planar limit, F^'^^ [T] is logarithm of the r-function of the dispersionless 
KdV hierarchy (i.e. Whitham hierarchy in the case of spherical Riemann surface). This is exactly the 
solution that was previously considered as relevant to 2d gravity [HJ |23l EH HOl E] . 

Throughout this section we denote 

9Tf 



U = ^ (4.1) 



and its dispersionless counterpart 

The Lax operator is 
The evolution is given by the flows 



92 ^(0) 



(4.2) 



L = g^d^ + 2U (4.3) 

^=.[LtL] (4.4) 

Here [•]4. denotes the differential part of the pseudo-differential operator (i.e. the " non- negative part" 
in the formal operator d). The appropriate solution of the Virasoro constraints is fixed by the string 
equation 

[L,M] = 2g (4.5) 

where 

N 

M = ^{2k + l)T2k+i L^"'/' (4.6) 

k=l 

4.1 An example: N = 2 case 

Let us assume that Z[T] = e^t'^l/^', which we "found" (up to some arbitrary function) in subsection 
12.21 solving the first Virasoro constraints, is, in addition, a r-function of the KdV (or the 2-reduced 
KP) hierarchy. The 3rd and the 5th equations of the hierarchy are 

3dT3~dTi\Y^^l2^ 



15 8X5 dTi\ 6 12 ''24 240 
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Picking up F from (j2.37p and inserting it into these equations leads to the Painleve-I equation on 
G{m) = F"im): 

g^G" = -6-G' + ^m (4.8) 

Solving it perturbatively w.r.t. g"^, one obtains series (I2.40p with the coefficients being recursively 
determineqj. The first two coefficients coincide with those in (|2.38p . From that equation it follows 

that F may also have an extra term C'r]2- However, one can check that actually C = using e.g. the 
equation 

= ^^l^ (4.9) 

where the residue is defined by 

oo 

res ^ ak{gdf = a_i (4.10) 

fc=— oo 

Thus, we have fixed all F^p^[T] completely, there is no more ambiguity and the curve and the 
densities are fixed. As we mentioned, F^'^^T] is such that the elliptic curve actually degenerates: 



y' = iBn) 1 (4.11) 



1 \/Ti-^TiTs+2T3 

i.e. the torus is pinched at the point y = 0, x = — ^ jr . Thus the curve is equivalent to 

..9 2 jTi-fTiTs-Ts 
the rational one Y = x — -^-^ tp . 

n 

A few first (multi)-densities are (here we consider the critical poinl|j case = 0, = ^ to make 
formulae more compact): 

viz) = -^(7t\ + z)y(z), Y\z) = z - 2Vt, (4.12) 



(012), , AVT,-z^-z, + ^Yiz,)Yiz,) 
P ^"^'"^^ = Y{zr)Y{z,){z,-z,Y = /'(N-)(^i'^2) 



(4.13) 

s=2v^ 



p^°^'^\zi, Z2, Z3, Z4) = ^- (Z1Z2Z3Z4 - 5 rj^^^(ziZ2Z3 + Z1Z2Z4 + Z1Z3Z4 + 2^2^:3^:4) + 

ATl'^Y^{z^)Y^{z2)Y^{z^)YHz^y 

+ 16Tl(ziZ2 + 2:12:3 + ZlZi + 22^3 + Z2Zi + Z3Z4) - 44rf^^(zi + 22 + 2:3 + z^) + 112 T^^) (4.15) 

naii)M = ^_ z-^VTi 

P 48. Ti Y^{z) ^ ' 



(211) r A 7 eiSTf - 503T^^^z + 204Tiz^ - UtI^^z^ + Az^ 

P^^ '{z) = — i i ^ 4.17 

9.210.2-7/2 yii(z) 



*Note that, strictly speaking, there are two solutions to (|4.8p corresponding to the two branches of in (|2.40p 
(or, equivalently, to the two branches of ^Ti - f TiTs in the expressions below) . But only one of these branches has a 
smooth limit T5 — > 0, i.e. permits the transition A'" = 2 — > A'^ = 1 considered in subsection 12.31 Therefore, only one of 
them is associated with the DV solution. 

^By critical point we mean the case when T2JV+1 — const and other Ti = 0, i > 1. Ti is as usual a variable. 



27 



Izfzl - Uy^{zjz2 + z^zi) + 35Ti{zl + z|) + 89TiZiZ2 - 182T'^^^{zi + za) + 284rf 



Thus, we checked that if Z = e^l^'^ satisfies the Virasoro constraints (for N = 2), then the following 
statements are equivalent: 

• Z is a T-function of the KdV-hierarchy 

• The elliptic curve y{z) degenerates into the rational one and the multi-densities p^Pl™) have no 
poles at the two marked points on the sphere which come from the double point singularity on the 
torus. 

Equivalently, one can say that the poles are only at y = points. Thus, the condition imposed 
on the curve determines F^'^^^T] and then the condition of canceling the singularities determines 

F(p)[T], p>0 (for general F^p^[T] the poles at x = '^'^s " exist). We wiU specify this 

statement for general case below. 

4.2 Quasiclassical limit of the KdV hierarchy as Whitham hierarchy 

In this subsection we review some features of the quasiclassical (or dispersionless) limit of the KdV 
hierarchy and its representation in terms of the generalized Whitham hierarchy |45t I46j . We work 
only quasiclassically and suppress the corresponding (0) superscript. In the quasiclassical limit of the 
KdV hierarchy, the momentum is just a commuting variable: gd P, because one can neglect all the 
commutators [P, ■] ~ g. Then, the Lax operator is just a function 

L = P'^ + 2u (4.19) 

which satisfy an additional constraint - the string equation: 

{L, M} = 2 (4.20) 

where {•, •} is the Poisson brackets {{P,Ti} = 1) and 



N N~l 

M = J2i'^k + l)r2fc+i = J2 MjP^^^' (4.21) 

k=l 1=0 



By [•]_!_ we denote the positive part of expansion at the vicinity of P = oo. and can be 
represented as follows (see proof in Appendix |A|) : 



Tk = -^resp=oo [PML-'^/^dP^ 
-resp=oo 



dF 

m 



(4.22) 



This defines our system as generalized Whitham hierarchy for the sphere parameterized by P. 

For generic N, the curve (|2.23p also degenerates to the rational one, when one imposes on Z an 
additional constraint to be the r-function of the KdV hierarchy. This is actually an obvious con- 
sequence of the equivalence of the curve (|2.23|) and the one appearing in the Whitham hierarchy, 
since the dispersionless limit of KdV (which describes F(''^[T] in (j2.23p ) corresponds to the Whitham 
hierarchy on the sphere. Hence, imposing on the solution of KdV hierarchy additional Vira- 
soro constraints quasiclassically is equivalent to total degeneration of the corresponding 
curve C2,2JV i given by (12.231) . This degenerated curve has the global parametrization 
z = L{pi y = M{P). 
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Indeed, let us denote 



N-l 



X = P2 + 2u = L; y=Y^ MkP 



2k+l 



k=0 



Then 



f{x) = (x - 2^) Yj ^k{x - 2uf 



, fc=0 



The curve C2 2N-1 defined earher reads exphcitly as 



(4.23) 



(4.24) 



y 



I / AT \ ^ I N-2 N 



z 



k~m—2 



{2k + l)T2k+i 



dF 



m=0 k=m+2 

To show that y^(x) = y^(x), one can use formulae (j4.22p 



\k=Q 



dt 



2m+l 



1 

'2k 



resp= 



dF 



1 

--resp=oo 



where dS = ydx. Then, expression (|4.25p can be combined into the sum 
1 



a+b>0 



resp= 



where 
Therefore 



/2 



resp= 



X 



1/2 



resT 



'-'/^dS 



dx 

y.a+1 



y\z) = [y-\z)]^=y\z) 



(4.25) 



(4.26) 
(4.27) 

(4.28) 



y{z) (4.29) 



(4.30) 



The inverse is also true. The condition of curve degeneration into sphere gives us A^ — 1 independent 
linear PDEs on F (pinching of all — 1 handles, or coinciding of the corresponding — 1 pairs of 
roots) which fix completely the ambiguity, the function F'^'^^ of — 1 variables. 

When all T^'s are fixed, one can completely determine all M^s and u. That is, u is defined by the 
equation 



N 

E 

fc=0 



(2fc + 1)!! 
k\ 



Tafc+i = 



(4.31) 



which has only A^ solutions. Thus, in generic case one obtaini 
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N-l 



l{x) = P{x)yr{x) P{x) = \/x - 2n, yr{x) = ^ Ms{x - 2i 



(4.32) 



'y(a;) of subsection 14. II is just P{x) 
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4.3 An alternative way to obtain multi-resolvents 

When Ti, . . . , T27V+1 times are turned on, there are solutions to the equations of the KdV hierarchy 
that satisfy the Virasoro constraints. The choice of solution corresponds to the choice of a root in 
equation (|4.31|) . — 1 of these solutions have a smooth limit when T2Ar+i — > 0, and one of them 
(corresponding to the root of (j4.3ip which goes to infinity) diverges. However, only one solution 
survives when all the times — > 0, A; > 3 reducing in this limit to the Gaussian model. Putting this 
differently, one may associate T^'s with r^'s of the Gaussian model in this case. Partition functions 
for different choices of T^'s are actually given by the same function in the sense that 

Z{T', t') = Z{T, t) = ZK{t) if T^fc+i + r^^+i = T^k+i + ra^+i = hk+i ■ (4.33) 

However, different solutions, i.e. different functions F[Ti, . . . ,T2Ar4.i], can be obtained from this one 
by analytic continuation with respect to variablefl] Ti, . . . , r27v+i- 

In this subsection we explain how to derive formulae analogous to (|3.34p and (j3.36p for multi- 
resolvents for the general DV solution considered in this section. This derivation is closer to the 
original way of getting ()3.34p in [20j (see also [l^). The crucial point for the derivation is to use 
specific moment variables [14J. The partition function in these variables can be easily obtained within 
the realization of the Kontsevich partition function as the highest weight of the Virasoro algebra given 
on the spectral curve (j3.6p . [6]. Within this approach, one makes a change of the local parameter on 
the spectral curvJ^ z — > G'^^^{z — 2u) (where G and u are some functions of times to be fixed yet) El 
that generates the change of times to the moment variables, 

1 nviz) + W{z))dz 

hm+i = 2;;r+r f -1 — (4.34) 

(2m + 1)G^ / (z - 2u)™+5 ^ ^ 



Then, the partition functions in old and new variables are related by the formula (see [6l sect. 3] for 
the detailed definitions and derivations): 

ZK{t) = G-Jie^^^^WZK(t) (4.35) 
This formula is correct for any functions G and u. We specify them so that 

n = rs = 0, (fsfc+i - ^Sk,i = t2fc+i) (4-36) 

Then, 

^(0) = UKK{t) = l<f <f P^'^^^\zi,Z2\2u)i(b{zi)+^Zl)){(P{z2) + <^>{z2))dzidZ2 

where p'^^^^^zi, Z2\2u) is given by formula ([3l8]) and (f){z) = ^rfcz'^+^Z^ $(z) = ^TfeZ^+^Z^ 2(f)'{z) = 
v{z), 2^'{z) = W{z) 

Introduce now, after [20j, the function 

,_ I (viz) + Wiz))dz _^ i2k + 1)11 , 



'^^We use it implicitly in what follows. E.g., formulae (|4.35|) - (|4.37p are known for the Gaussian branch. We derive from 
them explicit formulas for p'*''*^'s, e.g., (|4.52p . They are rational functions in u and Mk's (which in turn are polynomials 
in u and linear functions in T^'s). Once one have obtained these formulae for the specific solution to the KdV equations, 
one can make analytic continuation in Ti, . . . , r2jv+i to prove that they hold on over the branches. 

^^Note that the spectral curve considered in [6] is double covering of curve (|3.6p considered here. Hence, the local 
parameter ^ of [B] is related to the parameter z here as = z. 

^^As it will be seen later u — u|r=o coincides with u appeared earlier in section |4] 
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Then, from (fOij) it follows that 

S{u) = 

= -G 



da 

and, for k > 1, 



dS(a) (^-39) 



'^2fc+i = 2Fpr> (4.40) 

(2fc + l)!!(-5') — 



3 



where S^'^) = 5('=)(n). Now one can represent the operator V(z) in terms of the derivative with respect 
to u 

V{x) = f{x)n + ^{x) (4.41) 



2 



(x- 2u)" 



(4.42) 



where operator ^(a;) acts only on explicit dependence of S^^'^ on r's. Then, 

V(x).fW = /p(°|2)(x,z|2u)(</.(z) + cl>(z))dz - ^ L_|.^(0) (4_43) 

J (x — 2u) 2 ou 

The derivative of the two-point function splits 

|-p(0|2)(zi,z2|2u) = -3^^ (4.44) 

(zi -2u)2(z2 -2u)2 



and 

\J (x-2u)t J 



du" 



In this way one gets the well-known result, which is general for matrix models: the second derivative 
of the planar free energy depends only on ramification points (parameter u in our case) 

V(zi)V(z2).^^°^ = p(°l')(zi,Z2|2u) (4.46) 
Acting again with the operator V(z3), one gets 

V(zi)V(z2)V(z3).^(°) = ^— -T (4.47) 

'5'((zi-2u)(z2-2u)(z3-2u))2 

From the observation that [fi, V(x)] = with help of the relation 

k 

v(zi) • • • v{zk){s')-' = n\s')-' n f('k^ (^-^8) 

1=1 

one derives 

V(.i)-..V(z,).fW = -(-1|^)' 'ijllfiz.) (4.49) 

^ ^ i=l 

where, if one wants to get the expression for zero times r, one should substitute 



N 



k=s 



31 



A(2A; + 1)!! . 
u|,=o = ^^, u ^2fc+in^ = (4.51) 



k\ 

k=0 



which coincides with the string equation, (j4.3ip . Thus, we obtain the following generalization of (|3.34p : 

Let us note that = (2s + 3)Ms^i and M^, u have a nice description in terms of the curve (see e.g. 
g32D). From g3ZI)> and dOT]) one gets 

For higher genera it is probably too naive to expect expressions for multi-resolvents to be same simple 
as for genera and 1, but one can get (less simple) explicit expressions. If the variables u and 
Jk := ^C^) for /c> are considered as independent, then 



k=i 

d _ d ^ d 

k=l 



(4.54) 



with /W(x) = ^fix). Thus 



fc = l 

This operator is rather easy to apply to low genera, because in terms of the variables Jk, expressions 
for free energies are simple, for instance, 

1 / J4 29J2J3 7J| 



^(2) = _^ _ J^i + 1^ _ ^ (4.56) 
9216 V Jf Jf ^Jf. 



.4 79n T 5 ~'~ 1/1/1 7 6 ' 



192 Ji^ 720 Ji^ 144 Ji 

2 



1 Vn^) '^'•^^''^ 



ii52jf y "-^ 

For the Gaussian branch = M5q k, 2u = s, this gives 



/(4)(x) 105 



One can easily see that, in the generic DV case, /9^Pl"^^(xi, . . . ,Xm) can be expressed through the 
coefficients Mq, . . . ,MAr_i and P{x) only. From formula (j4.55p it follows that the singularities in Xi 
can only be of the form f^''\xi)\r=o, i.e. l/P^'^+'^(xj). 

Here we list several first densities (for general N) (using variables n, Mq, . . . , M^v-i instead of 
,r2Ar+i): 

p(o|i)(x) = -W{x)/^/^ - P{x) ■ yr{x) (4.59) 
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^(0|2) 



(xi,X2 



1 



P 



(°I')(X1,X2,X3 



X2) 



P{XI)P{X2) 
1 



^/XlX2 



MoP^Xi)P^X2)P\xs) 



1 



+ 



1 



1 



+ 



1 



-1 Mo-MiP2(a;) 



Ml 
Mo 



8Mo2 



P5(x) 



(4.60) 



(4.61) 



(4.62) 



(4.63) 



Note that the dependence on x and u enters only through the difference x — 2u = P'^{x) (analogously 
to what we had in the = 1 case (see formula ()3.8p ). again except for the non-meromorphic parts in 
yo'-'^l^^ and /9^'^|2)), This type of solutions to the Virasoro constraints is sometimes referred to as one-cut 
solutions (for the obvious reason). 

4.4 KdV/DV solution and absence of poles 

Taking into account what was said at the end of the previous subsection, one can state that if Z{T) 
satisfies the reduced Virasoro constraints Lq? then Z\T] = Tkp <^=^ the curve C2^2N~i 
degenerates to a sphere and for all p, m have no poles at marked points on the 

sphere which came from degenerated handles, i.e. at the zeroes of yr{x). 

Proof of the genus-zero (regarding degeneration of the curve) part of this conjecture is presented in 
s l4.2l It is actually enough to check the condition for higher p only for the one-point functions p^'P^^\x), 
because when performing iterations (j2.20p the singularities just can not arise in the multi-densities 



P 



{p\m) 



(xi,..,Xm) due to the symmetry between xi,..,Xm- Then, for each p we have the following 
condition: the numerator of p^P^^\x) obtained through (j2.20p . which is a polynomial in x, should be 
divisible by yr{x) in the denominator. In the general case, the remainder is a polynomial of degree 
N — 2, and thus one has N — 1 equations, which are the first order linear PDEs on F^p\ Hence, the 
function F of — 1 variables is completely determined. 

The absence of the poles in all is equivalent to vanishing of the integrals 



p^^\x)dx = 0, 



Vz = 1 . . . - 1 



(4.64) 



where the ^-cycles encircle handles (which are actually pinched): 





no poles here, 
cycles can be 
contracted 



and pl^^ = YlpLo g'^^ P^'^^^^ is the total one-point function. Therefore, one can actually forget about this 
points and the initial curve, and work in terms of the reduced curve P{x) = \/x — 2u. 
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4.5 A particular case: Chebyshev spectral curves 

There is a special choice of values of which is usually referred to as conformal background 
This choice leads to the Chebyshev curveJ^. Indeed, if values of times T/. = Tk{X) are given by 

f;T,.,,e'=+V.(^^ = '-^i-Xf-'^' [K,.y,i-X/mU (4-65) 
then there exists a solution u = —j of the string equation ()4.3ip such that 

x = P'-^A = ^T2(P/VA) (4.66) 

y = A^-V2T2Ar_i(P/VA) (4.67) 

or, equivalently 



X [j ) 



(4.68) 



We denote here by T=k the Chebyshev polynomials. The choice of times (I4.65P defines the conformal 
background. 

In order to prove this claim, one can set A = 1 without loss of generality. According to what is 
said at the end of s l4.21 the space of curves of type (|4.23p - (|4.24p covers N times the space of T/^'s. 
Then, it is enough to prove that the times given by the Chebyshev curve y = Tjv-i/2(2x) are such 
that (I4.65P satisfies. Taking into account that T2k+i = — 2k+i ^^^x=oo [x~''~^^'^y{x)dx^, it is correct if 



N 



^ (T,)fcr(fc+l/2)x"-i/2 = (_i)-^.i^^(_l/^) (^,^Ar_i/2) (4.69) 



fc=— oo 



is correct. We denote here by K^: the modified Bessel functions of the second kind and Ti/(x) = 
X^fc(T!/)fex'^'~^/'^. Using the integral representation for the T-function and deforming the contour, one 
can rewrite this equation as 



+00 

T^{T)e-^HdT = {-iyi^K„{-l) {l = l/x>0) (4.70) 



Then, integrating by parts in the l.h.s. and using the integral representation for the Bessel function 
in the r.h.s., one comes to 

+00 -If 

T;(r)e-^' dT = - e-2(*+7) i^t^-^dt (4.71) 



1 2 jc* 

where contour C goes from 0+ to +00 encircling one time. One can check this equality just by 
changing the variables t + l/t = t , t = t± = t ± a/t^ — 1 and using T^(r) = (t+ + t!l)/2. 
Note that in the general case if one has a curve 

Tp(x) = Tg{y) (4.72) 

then due to the composition property of the Chebyshev polynomials, Tp o Tg = Tpg the curve can be 
parameterized in the obvious way 

X = T,(P) y = Tp(P) (4.73) 



In [49] it was shown how the Chebyshev curves emerge in Liouville theory. 
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5 Generalized Kontsevich Model and duality 



In this section we review some basic aspects of GKM. The next section contains some more detailed 
analysis of the p = 3 case and includes finding the curve, the first densities and action of the p — q 
duality on them. 

5.1 Generalized Kontsevich Model 

Partition function of GKM is given by the following matrix integral ([10]): 



It is a function of the Miwa variables 

Tk + Tk = tk = ^<5p+i,fc + ^TVA-^ (5.2) 



Actually it depends only on times t^, k ^ O(mod|)). The Ward identities in this case are the W- 
constraints fpkllTn]) 

Wi^^Z = 0, k = 2..p, -k + l (5.3) 

(2) 

For k = 2 these are the Virasoro constraints WA = Cn- One can define the densities analogously to 
the p = 2 case 



p(^l™)(zi, . . . , z„) = V(zi) • • • V(z2)^(^)|.=o, V(z) = (5.4) 



ifc=0 

In particular, 

p+<i ^ 

p(o|i)(x) = -VF(x) + y(a;), W{x) = Ykxv~^Tk (5.5) 

k=l 

and one expects y{x) to satisfy an algebraic equation 

Pp,,(y,x) = (5.6) 

defining the spectral curve. The polynomial Pp^q{y,x) is of degree p in y and of degree q in x. The 
GKM is known to be relevant for describing the {p, g)-model of 2d gravity when the first p + q times 
are turned on (T„ = 0, n > p + q). The well known solution to ()5.3p is given by the r-function of 
the p-reduced KP-hierarchy ([44^ I23j). The spectral curve in this case degenerates to the rational one. 
When one further specifies to the conformal background the curve has the form of 



Tp{y)=T,{x). (5.7) 



5.2 p-q duality 

The p — q duality is obvious when one formulates theory in terms of the Douglas equation [H] gener- 
ahzing (j4.5p 

[L,M] = 2g (5.8) 

where L and M are differential operators of orders p and q respectively. However, when one formulates 
the theory in terms of p-reduced KP-hierarchy (where L plays role of the Lax operator) or in terms 
of the matrix model, the duality becomes implicit. In an explicit change of time variables T T 
of the hierarchy which connects {p,q) and {q,p) models was found. In [T9j the p — q duality was 
considered from the point of view of the GKM. 

Within our approach, one first of all expects that the duality relates the functions F'P'''' and F^'i'P) 
(generalizing those considered in s ]2.2p which describe ambiguities in the solutions of the W-constrains 
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(so that fixing them is equivalent to fixing the solutioio). They both depend on a certain number 
(equal to the genus of the spectral curve) of combinations of times T, consequently the change of times 
sliould map the set of these combinations for the {p, q) model to that of the {q,p) model. Then (after 
identification of F'''"'*' and F'^'P') one expects that the curves for the two models should coincide, that 
is, there is an isomorphism C(p,q) C(q,p): 

Pp,g(y<P'''', x'P''^') ~ Pg,picl)*y^'''^\ c^*x^^-^^) (5.9) 

In simple cases, this is just the interchange x ^ y: 0*?/''''^' = j;'^'''', (p*x'-^'^^ = y'p-'i'. Note that, as soon 
as we talk about solutions to the W-constrains only, without further specifying them as KP solutions, 
in variance with [n\ \TE[ [T9] we deal with a more general, "off-KP" duality. 

6 Generalized Kontsevich matrix model (p = 3) 
6.1 Some general formulae and definitions 

Most of what comes up below is quite analogous to the p = 2 case, hence, we present it very briefly 
(one can find some parallel consideration in |50j). 

>vi^'-^»4fe ^ '''''' ^ ^ E 

\ i+j=~3n i—j=—3n i+j=3n I 

\ p+q+r=~'in p+q~r=~3n 

,2 V- ,_^,9^ V 

+5 PPdtgdtr 3 dtpdtgdtr 

p—q—r=—3n ^ p+q+r=6n ^ ^ 

CnZ = 0, n ^ -1 (6.2) 

W!^^Z = 0, (6.3) 

= 0, = (mod 3) (6.4) 
The shift of times: tk = Tk + Tk, / for /c = 1 . . . ((7 + 3). Thus, it can be referred to as (3,q) model. 



oo 



V^(^)-E^a;^ ^ = 1'2 (6.5) 



(x) = ^(3A; + i)xSfe+i Wi{x)= {^k + i)x^T3k+i i = l,2 (6-6) 

fc=0 fc=0 

Pi{x) ■.= V,{x)J' i = l,2 (6.7) 



The loop equations are 

°° 1 2 2 

E -^^n^ ^ = 2(^1 + ^l)(^2 + ^2) + 1^ + Plix)p2ix) + y (Vi(x)/52(X) + V2(x)pi(x)) + 



n=-l 



+P- {(ni(x) + Wi{x))pi{x) + (n2(x) + W2{x))p2ix)} = 



'Still, fixing the solution is equivalent to fixing periods of all the multi-densities. 
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n=-2 ^ 



^^^Ix-lnez+l : (s'V -iu + W))^ :} 



fc=0 



3A;+i 



V(x) = Vi(x)x"^/^ + V2(x)x~^/^ = /3i(x)x~^/^ + p2(x)x"2/^ = V(x)J^ 

W{x) = VFi(x)x-2/3 + W2(x)x-^/3 



F[Ti, . . . , Ts+g] — JF|t-=o 



p(P\^\z,, ...,Zm)= Viz,) • • • V(z2)^(^)|r=0 



6.2 Some first densities 

Let 

(o|i) 
Pi 



where 



From the definition of p^^^^^ one has 



W2 + y2X^/^ 4°'^^ = -VFi + yix2/3 
p(0|l)(x) = +y(x) 

2/1 + y2 = y 



(6.8) 
(6.9) 
(6.10) 
(6.11) 
(6.12) 



(6.13) 
(6.14) 
(6.15) 



{x-'=-^/3y(a;)dx} 



-res,=oo {x'=+*/='y(x)dx} (6.16) 



r=0 



Substitutin£2| ([613]) into the loop equations for genus zero one obtain the following equation for 
y defining the C(p^3_q) curve: 



/ + ?>A{x)y + B{x) = 



where A and B are the following polynomials: 



Aix) = -P+ {Tyi(x)W2(x) + [Tyi(x)(Vi(x)F(°)) + W2(x)(V2(x)fW)] } /, 



B{x) = -P+ ^Wf/x + Wi + 3 P^i(V2F(°))Vx + Ty2(ViF(°)) 



+ 



+3 



VFf (V2F(0))/x + W|(ViFW)] } / 



^^What follows is actually true only for q < 7. 



(6.17) 
(6.18) 
(6.19) 
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There is also the following identity: 



yiy2 = -A (6.20) 

which, together with yi + 2/2 = V-, allows one to express yi through y and A. 
Some formulae which are useful for computing the higher densities are 

V^{x)p-{ui{z)h{z)] = -3xl5,x-| I Kx)^ - Kz) I ^ p-{^.(^)v,(x)/.(z)} (6.21) 



I B' + 3A' y 

3 y'^ + A 

An explicit expression for in a certain special case can be found in s l6.6 



6.3 Solving reduced constraints and (3,2)^(2,3) duality 

Now c 
ingly: 



(3) 

Now consider reduced Virasoro constraints for q = 2 coming from WI21 £o-constraints accord- 



lr| + 2TlT, + lOr^Ts^ + 8T||| + |r|g = (6.23) 

dF OF 

2TiT2 + 4T4— + 5T5— = (6.24) 
oTi dT2 

o2 dF OF OF dF 

— + ri h2To h4r4 hST^ = (6.25) 



The general solution is as follows 



^ ^ ^ , g2 1 TiT2^ 4 T2^r4 2 T^Ti^ 16 T/Tir2 32 74^2 

F[T,,T2.T,,n] - "Y^logTs - - -J- + - ^ - - + — - — 

128 T^^Ti 512 r4'^r2 4096 T^^^ ~ ( 125 TiT^^ - 100X2X4X5"^ + 16 T4'^\ 
' + F 



9375 46875 Tg^ 5859375 Tg^ ^ 54 . 34/5 . 2-5/5 . (_r5){i6/5) J 

where F is an arbitrary function. This is the (3,2) model. The (2,3) model corresponds to the A'"'*' = = 
2 case considered in subsections 12.21 and 12.31 Let us denote Ti = T^^^^^ , T3 = Tg'*^^' , T5 = T^^^^^ . 
Using the transformation^^ (the same as in jl7]): 

T,=T,-1^-A^ + 1^^-1^ (6.26) 



5 Tg 20 Ts 25 Tg^ 25 T; 



5 

A rp 2 

fs = T,--'^ (6.27) 
5 Ts 

r5 = -^r5 (6.28) 



one can see that 



TiTs - ^T^ 125 Ti Ts^ - 100 T2 Tg^ + 16 T/' 



^0 ^^7^ ^7(167^^ (6-29) 



■^5 



^'^The variable T3 can be considered as auxiliary. All (3,2) quantities are actually independent on it. For instance, one 
can choose ^3 = that the last two terms in the r.h.s. of (|6.30p cancel each other. 
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and, thus, the arbitrary functions in the two models can be identified: -F{2,3) = -^(3,2)- Moreover, one 
gets the same relation as in jl7] : 

u,^ATi.T2,T,,n] = «(2.3)[Ti,r3,r5] + i ^ - ^ J (6.30) 

_ Q^-^(3,2) _ d'^F{2,3) /«Q1\ 

■"(3,2) - Qrp2 ""(a^S) - Qj,2 (^D.Ciij 

Note that we used so far the W constraints only and did not assume that Z is a r-function of the 
KP-hierarchy. Therefore, this is an off-KP duality. 

6.4 Duality between (2,3) and (3,2) curves 

For A^'p"-*' = 1 (i.e. the (3, 2)-model) we have the curve 

+ 3A{x)y + B{x) = (6.32) 

A{x) = -2OT4T5 ■ X - (STiTs + 8T2r4) (6.33) 



B{x) = -125Ti ■ - (ISOTaTl + 64T|) • x - {GOT^n + ^STiT^ + TSTg^— |^) (6-34) 



dTi 

This is an elliptic curve. For A^(p=^' = 2 (i.e. the (2, 3)-model) we had the following elliptic curve: 

_ dF^^^ 

y2 - 25r| • x^ - 3OT3T5 ■ x^ - (lOTiTg + 9r|) • x - {GTiTs + lOTs— |^) = (6-35) 

After the change of variables T ^ T and the identification -F(2,3) = -^(3,2) described in s l6.3^ one 
can check that the j-invariants of the curves coincide with each other: 

/2>3) = ^-{3,2) (g_36) 

i.e. the curves are isomorphic. At the critical points, the isomorphism is given just by 

3.(2,3) _ y(3,2) y(2,3) ^ ^(3,2) (g_37) 

Further details on equivalence of these Riemann surfaces can be found in Appendix B. 
Note that this is again the off-KP duality. 

6.5 Duality between p*^°l^^'s 

The isomorphism between the elliptic curves C(3^2) ^(2,3) is given at the level of coordinate functions 
by the following linear (!) explicit expressions 

</,*y(-3) = • x'^-^) + • + 2 + (6.38) 

<^V^,3)^+1 1 .y(3.^) + |^ (6.39) 
5 Ts 5 r5 

Using this, one can see that 

mod exact form) (6.40) 



„(0|1) 

A^(2,3) 



(OH) 

"(3,2) 



where we regard 

p(0|l) = ^(o|i) 

{x)dx as a 1-form, and [p^'^'"'^^] j^^^,^ is its meromorphic part (without the 
W^(x)-term, which simply cancels the singular part of the expansion of y w.r.t. x nearby x = co and 



39 



thus, given the local coordinate x, can be recovered easily). Therefore, the meromorphic differentials 
y (3, 2) ^^{3,2) y{2,3)^^{2.,3) jjffgp ^ meromorphic differential with zero periods. This can be also 
expressed as 



„{0|1) 

"(3,2) 



JOID 

"(2,3) 



JOID 

"(3,2) 



"(2,3) 



(6.41) 



A, 



Bi 



where j4j, Bi are the cycles on the C(3^2) (we assume that the contours do not encircle the ramification 
points X = and x = oo produced by VF(a;)-terms in in both cases). 



6.6 On in (2,3) and (3,2) models 

The isomorphism (j), of course, connects the holomorphic differentials on the curves: 

For the sake of simplicity, in what follows we work at the critical point (only Ti and are switched 
on). Then one has 

dF 



oil 



y 



(2,3) 



-2 • 5 • T. 



^^(3,2) 



(y(3,2))^+^ 



In the p = 2 case, there is a decomposition of p^^^"^^ 

,(0|2) 



- .(0|2) , J0|2) _ (0|2) 
~ Phol ^ Pglob Ploc 



(6.43) 

(6.44) 
(6.45) 



(0|2) 
P(3,2),hol 



dxdz 



^ STf (y2(^)+^)(y2(^)+^) 



(6.46) 



(0|2) 
P(3,2),glob 



dxdz 



4(x - zy{y'^{x) + A)(y2(2;) + A) 



[2 • 3 • 52 . TfTl + 2 • 3 • 5 • T1T5 • y(x)y(z) + 



-y{x) { 2WTi^^ + 5'Tiiz' + 2x.)} - y{z) {2's'5'Ti^^ + ^'tUx' + 2xz) 



Note that pf}^) giob^ P%K) hoV Pt-itjioc have the same defining properties as p^^\^^ discussed in s J2.4l One 
can also try to formulate it in terms of the correlators (see s l2.5p . however, here this is more difficult 

^(3,2) 1 

to find an appropriate local CFT. The scheme of cuts of the covering C(3,2) — 



J0|2) 



(0|2) 
P (3,2), loc 

J0|2) 



[x^/^ {2z + x) + 2^/3 (2x + z)] dx dz 

4^2/3 ^2/3 (^^ _ _2)2 



(6.47) 



(6.48) 



12) 



IS 
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Here [a] and [b] are the standard homology basis on the torus. If one considers, say, the coUection 
of fields (Xi, X2, X^) on the sheets, there is no linear combination of them which diagonalizes the 
monodromies around all points simultaneously. 
For the (2,3) model one has 



{0|2) 
P(2,3),hol 



d'^F dxdz 



(0|2) 
P (2,3), glob 



The singularity at x = z: 



dxdz 



{x - zfy{x)y[z) 



.{0|2) 



22 5 y^^^y^^^ 

1 BF 

- {5'Tixz + lOfiTs} ix + z) + 5f5^ 

fx + z] dx dz 



<.2,3),loc 4^1/2 ^1/2 (^_^)2 



(0|2) 
P(3,2),glob 



(0|2) 
P(2,3),glob 



P, 



P, 



i|2) 3 dxdz 

(3.2) ,loc ~ 2(x-z)2 

(o|2) 1 dxdz 

(2.3) Joe ~ 2(X-Z)2 



Both p^^^^3) and /0(3 '2^ can be constructed in terms of the coverings C^2,3) 
The curves are isomorphic C(3,2) — *■ ^-(2,3) but, of course, vr(3 2) 7^ '^{2,3)4'- 



.(0|2) 



T(2,3) 



(6.49) 
(6.50) 
(6.51) 

(6.52) 
(6.53) 



»^ andC 



T(3,2) 



(3,2) 
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A Derivation of 022]) 

Lemma A.l There exist the following formulae for Tj^ and 



--resp=oo 
k 



resp=^[PML'''^dP] 



(1.54) 
(1.55) 



Proof : Let us denote 

^:=f;P^r,.«„' (1.56) 

Then, one can rewrite (|4.20p as = 0. From (j4.2p and the dispersionless counterpart of (j4.4p it 
follows that 
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oT2k+i k\ dTi 

Then 

thus H = (note that we have also checked in the way that this equation respects all the KdV- flows). 
One can see explicitly that the equation H = can be represented as follows 

Ti = -res |PML-i/2^p| (1.60) 

For k > the equations 

T2k+i = -^^res{PML-(2^-+i)/2dp} (1.61) 
follow straightforwardly from formula (j4.2ip for M. Thus, it remains to verify that 



§^ = -.e. {pM.-./=.p} = X:(-irM._..-^^- P'-;f'^P-"" (1.62) 

oT2k+i ^ '7^1 {s + k + l)\ 

Differentiating it w.r.t. Ti (note that, doing this, we lose nothing, because in the KdV hierarchy only 
the quantities containing, at least, first derivatives of F in Ti make any sense) and taking into account 
that, from the string equation ()4.20p . one obtains 

^ = (2. + l)M.|!i , Mo|^ = -l (1.63) 

one finally arrives at 

_ (2^ + 1)!!^,+! 



dT2k+idTi {k + l)\ 

which is the same as (jl.57p . □ 

For reference, we also write down explicit formulae for T^, Ms and ^§^- 

1 

T2k 



(1.64) 
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B On equivalence of Riemann surfaces 
B.l y'^ — x'^ — 1 and other hyperelliptic curves 

The complex curve = — 1 is an ordinary Riemann sphere with coordinate z, where x = + 1 /z) 
and y = ^(z — 1/z). There are no holomorphic differentials on it, and the kernel, the bi-differential 
with a double pole on the diagonal can be represented in a variety of ways: 

s dzidz2 1 dxidx2 ( xiX2-l\ 1 dyidy2 ( yiy2 + l\ 
[zi- Z2y 2 {xi - X2y V 2/12/2 / 2 (yi - ys)^ V ^1X2 J 

For the generic hyperelliptic curve, = Pn{x) = anx"^ + . . . + oq 

„ 1 dx-[dx2 I P„(xi,X2)\ 

^=o ( \2 h+ 2.69 

2 (xi - X2Y \ ym J 

where the new polynomial P„ is defined by the three conditions: Pn{x,x) = Pn{x), symmetricity 
Pn{xi,X2) = Pn{x2,xi) and restricted growth at infinity, Pn{xi,X2) ~ x^^'^^^^'^K This means that for 
even n = 2m 



,,'"-1 I ni 
P2m{Xl,X2) = a2mXi X2 + a2m-l ^ V 

™m™m— 2 I m— 2 m ^ 



(2.70) 



(XX H~ X X \ 
a2m-2^-^ ^-^ ^ + (1 - a2m-2)x'^'^x'^-^j +... 

while for odd n = 2m — 1 



„m„m— 1 I m— 1 m 
P2m-l{Xl,X2) =a2m-l 7i 



Extra a-parameters describe possible bilinear combinations of holomorphic differentials which 

can be added to B, for example, a!2m-2 is actually a coefficient in front of x^~^x^~^(xi — 0:2)^. This 
ambiguity can be fixed by specifying ^d-periods of the kernel B. In particular, the Bergmann kernel 
corresponds to the trivial A-periods 

B = (2.71) 



/ 



B.2 

This curve is an ordinary torus with extended discrete Z2 x Z3 symmetry. In the dual coordinates 
X = y, Y = X it acquires the usual hyperelliptic form = + 1, and the (un-normalized) 
holomorphic differential and the Bergmann kernel are 

In coordinates x, y one has instead a representation as a triple covering with the three ramification 
points at a; = ±1 and x = 00. The local coordinates in the vicinities of these points are 



Accordingly, the holomorphic differential and the Bergmann kernel acquire the form 



(2.73) 



2dx ^ _ dxidx2 / 2/1 + 2/12/2 + 2/2 V 1 (2/1 + y2)dxidx2 ^274) 



3 y2 ' (^xi - X2y V 3yiy2 / 9 yjyj 
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In the first term in B the numerator vanishes when 1/2 = ibe^'^*/^7/i and cancels the unwanted poles at 
the points xi = X2 with yi ^ y2- The second term serves to cancel poles at infinity. 
The modified Begrmann kernel, with the pole at X2 = —xi is 



^ 1 dXidX2 f XfX2 + X^Xl + 2 



2 (Xi - X2Y \ 2Y1Y2 



dxidx2 (yi+ 2/12/2 + yi ^ ^ _^ 1 (yi + y2)dxidx2 



(2.75) 



(xi + X2Y V 32/12/2 J ' 9 ylyl 



B.3 r = Mx)y + C2{,x) from = p^^x) 

Here ^i(x) = aux + oio, C2(x) = C22X^ + C21X + C20 and P2.{X) = p^^X^ p^^x'^ ^ p^^x + psQ. 

The equation y^ = Ai{x)y + C2{x) is quadratic in x, and the function x{y) has the four order-two 
ramification points at 2/ = 00 and at the three roots of discriminant in the expression 



-(c2i + any) ± V(c2i + aiiyY - 4c22(c20 + Qio2/ - 2/^ 

2C22 



(2.76) 



At infinity 



and 



Therefore, 



(2.77) 



rr-l, 2/-^ (2.78) 



Y = ux + vy + X = py + q (2.79) 



The holomorphic differential ^ ~ ^^^'^^^^ does not have poles at zeroes of the denominator, provided 
these zeroes are located exactly at the ramification points, where y has double zeroes. This means 
that actually 

Y = ux + vy + w '-^ 2c22X + a\\y + C21 (2.80) 
Substituting this expression into Y^ = P3{X) one obtains 

{P33P^) y^ + (3p33P^g + P32P^ - all) + {^P33PQ^ + '^P32Pq + P3lP " 2aiiC2l) 2/ - 4oiiC22X2/ = 

= 4C22X^ + 4C22C21X + Cgi - P3{q) 

which coincides with y^ = Ai{x)y + 6*2(2;) provided 

2af,=p2p^'(g) 

P33P^ = 4C22 
Psiq) = (?2\ - 4C20C22 

2aiiC2i - 4aioC22 = P-P3(q') 

These equations can be used to define p^q and ^i(x) for given P3(x). The remaining freedom is 
X ^ ax + /? and it can be used to fix two coefficients in C2(x). 
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